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Abstract. We present our results on neutral pion production in pp and dAu 
collisions at RHIC energy. Pion spectra are calculated in a next-to-leading or- 
der (NLO) perturbative QCD-based model. The model includes the transverse 
component of the initial parton distribution ("intrinsic kr")- We compare our 
results to the available experimental data from RHIC, and fit the data with high 
precision. The calculation tuned this way is repeated for the dAu collision, and 
used to investigate the interplay of shadowing and multiple scattering at RHIC. 
The centrality dependence of the nuclear modification factor shows a measurable 
difference between different shadowing parameterizations. 

1 Prologue 

With this paper we would like to pay tribute to Jozsef Zimanyi (Jozso) who passed away in 
September 2006. Jozso was an enthusiastic supporter of the Hungarian participation in the 
relativistic heavy ion program at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory. He had a detailed vision about quark-gluon plasma formation at RHIC, 
and he was eager to confront the real world of experimental data with his ideas. In the last 
decade of the construction phase of RHIC he predicted many observables, including hadron 
numbers and ratios [1]. He realized that the era of high-p^ physics and QCD is coming to 
nuclear physics, because the RHIC detectors extended the measurable transverse momentum 
region much beyond the SPS values of pr — 3 — 4 GeV/c. He followed the progress in this 
subfield and was always happy to discuss related questions. 

Jozso joined the Hungarian experimental PHENIX group from its foundation and supported 
continuously the Hungarian activity by full heart and by his authority. The first PHENIX 
paper on which his name appears is entitled "Absence of suppression in particle production at 
large transverse momentum in ^/snn— 200 GeV in d + Au collisions" [2], and has become a 
fundamental reference in the quest for understanding hadron suppression and jet energy loss in 
deconfmcd matter in AuAu collisions. Jozso is among the authors of many other dAu papers, 
including a detailed analysis of centrality dependence of the nuclear modification factor in dAu 
collisions [3], the results of which will be used in our present theoretical analysis. The results of 
Ref. [4] about identified particle production and Ref. [5] about charge particle production will 
also be commented on. 

Our group has performed theoretical investigations connected to the pp, dAu, and AuAu 
collisions at RHIC energies in the framework of a RMKI-ELTE-Kent-Columbia collaboration. 
We would like to posthumously thank Jozso with this paper on dAu collisions for his continuous 
interest in our work and for his support to keep our theoretical collaboration vital and successful 
in high-energy nuclear collision research. 
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2 Introduction 

Recent experimental data on mid-rapidity tt° production in the transverse momentum region 
2 GeV/c < pt < 15 GeV/c measured at the Relativistic Heavy Ion Collider (RHIC) in AuAu 
collisions at y/s = 130 AGeV and 200 AGcV display surprising evidence of new phenomena 
appearing in hot dense matter. Here we focus on the property that the n° transverse momen- 
tum spectrum shows a strong suppression [6-8] in central collisions compared to theoretical 
expectations based on binary collision dynamics in perturbative quantum chromodynamics 
(pQCD) calculations (see e.g. Refs. [9-11]). A similar suppression pattern was observed for 
charged hadrons [12]. Introducing final state interactions, especially jet energy loss [13-19], this 
suppression pattern can be reproduced both at 130 AGeV [20,21] and 200 AGeV [22,23]. The 
entanglement of initial and final state interactions in AuAu collisions and the complexity of 
the theoretical calculation generated an increased interest in simpler reactions, such as the dAu 
collision. In this case final state interactions play a minor role, and the initial state interactions 
(initial multiscattering, shadowing) can be investigated cleanly. The obtained dAu data from 
RHIC [24-26] stimulated further interest in theoretical calculations of pion production in this 
reaction. 

In a pQCD-based leading order (LO) parton model a systematic analysis of existing pA 
data has been completed [11] and our calculation will be based on results obtained there. A LO 
calculation of dAu was reported in Rcf. [27]. Here we perform a dAu computation in next-to- 
leading order (NLO), utilizing recent developments for pp collisions [28]. Successful calculations 
based on the physics of gluon saturation have also been carried out for this system [29,30]. 

In this paper we take advantage of the availability of detailed (impact-parameter selected) 
data and calculate the impact-parameter dependence of nuclear modification factor RdAu in 
our pQCD based model. 

The paper is organized as follows. In Section 2 we review the NLO pQCD-improved parton 
model augmented with the intrinsic transverse momentum distribution for pp collisions. We 
compare calculations to pp results at midrapidity at y/s — 200 GeV and use these data to tune 
the parameters of the model. In Section 3, we discuss calculational details for the d + Au — > 
7r° + X reaction at y = in the transverse momentum region px > 2 GeV/c. In Section 4 
the RdAu nuclear modification factor is extracted from the data, and we investigate its impact 
parameter dependence. In particular, we study simple shadowing prescriptions without and 
with impact parameter dependence. In Section 5, we discuss our results. 



3 The pQCD improved parton model with intrinsic k T 



The invariant cross section for neutral pion production in a pp collision can be described in 
the NLO pQCD-improved parton model on the basis of the factorization theorem as a convo- 
lution [28,31,32]: 
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where we use a product approximation for the parton distribution functions (PDFs), 

f(x,k T ,Q 2 ) = f(x,Q 2 )g(k T ) . 



(1) 



(2) 



Here, the function f(x, Q 2 ) represents the standard NLO PDF as a function of momentum 
fraction x at factorization scale Q, da/dv is the Born cross section of the partonic subprocess 
ab — > cd, K ab c (s,v,w,Q,Q R ,Q F ) is the corresponding higher order correction term, and the 
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fragmentation function (FF), D^(z c , Q%), gives the probability for parton c to fragment into a 
pion with momentum fraction z c at fragmentation scale Qf- We use the conventional proton 
level (S, V, W) and parton level (s, v, w) kinematical variables of NLO calculations (see details 
in Refs. [28,31,32]). In our present study we consider fixed scales: the factorization and the 
renormalization scales are connected to the momentum of the intermediate jet, Q = Qr = K-p q 
(where p q = pt/z c ), while the fragmentation scale is connected to the final hadron momentum, 
Qf = k ■ pt- The value of k can be varied in a wide range, k € [0.3, 3]. 

Our NLO calculation includes the initial transverse- momentum distribution .g(k/r) of partons 
("intrinsic kx"), along the lines of Refs. [9,11,28,33]. We demonstrated the success of such a 
treatment at LO level in Ref. [11]. In our phcnomcnological approach the transverse- momentum 
distribution is described by a Gaussian, 

fl (k T ) = ^e"*/<*> . (3) 

Here, (k T ) is the 2-dimensional width of the kx distribution and it is related to the magnitude 
of the average transverse momentum of a parton as (k T ) = 4{kx) 2 /TT. 

In our LO and NLO investigations we use GRV LO PDFs [34] and the MRST(cg) NLO 
PDFs [35] for parton distribution functions. For fragmentation functions we use the KKP 
parameterization [36], which has LO and NLO versions. An advantage of the GRV and the 
MRST sets is that they can be used down to very small scales (Q 2 w 0.25 GeV 2 /c 2 ). Thus, 
they provide reasonable calculations at relatively small transverse momenta, pr > 2 GcV/c at 
our fix scales. 

Figure 1 displays our LO and NLO pQCD results at y/s = 200 GeV for pp -> ir°X [37] at 
different scales. In the LO case (dashed line) we used k = 2/3 for both the factorization and 
the fragmentation scales to reproduce the data at high pt within error bar. (This scale value 
is higher than k = 1/2 used in Ref. [11] due to a readjustment triggered by the availability of 
precise data at high pt on n° production in place of the earlier UA1 data on charged hadrons.) 
To minimize the difference between RHIC data and the model at low px a (k T ) = 2.5 GeV 2 /c 2 
is included. We obtained quite good agreement in LO (see open squares in the data/pQCD 
ratio, lower panel). 

In the NLO calculation we use the scales Q = Qr = (4/3)p g and Qf = (4/3)pr to reproduce 
the central value of the experimental data at high px (dash-dotted line in the upper panel). 
However, in the window 2 GcV/c < px < 4 GeV/c this choice underestimates the data by a 
factor of 2 (see filled triangles in the lower panel). These calculational results are in agreement 
with the calculation reported in Ref. [37], albeit with a different set of scales [38]. Since we 
plan to investigate multiscattering and the Cronin effect in dAu collisions, which are most 
pronounced in the above momentum window, we consider an NLO parameter set with scales 
Q = Qr = (4/3)p 9 and Qp = (4/3)pt, together with an intrinsic kr value (k T ) = 2.5 GeV 2 /c 2 . 

Fig. 1 indicates that we obtain good agreement both at high and low pt with this choice (full 
line in the upper panel, filled squares in the lower panel) . One can see that the results with this 
NLO parameter set and our LO results are close to each other, even though the scales differ by 
a factor 2. This is due to having the same width of the intrinsic kx distribution. The physically 
realistic nature of the value found for the width is corroborated by results from di-jet production 
at lower energies: jet-jet correlations yield a similar value for (k T ) at ISR energies [39]. While 
recent measurements of jet-jet correlations in pp collision at RHIC energies [40,41] may further 
clarify the properties of the transverse component of the PDFs, these paramctcrizations fixed 
in pp collisions will provide a solid basis to investigate nuclear collisions. 

The precision of the NLO calculation at lower px could be improved introducing non- 
Gaussian or pT-dependent intrinsic kx (see e.g. Refs. [42,43]). In the present analysis we keep 
a constant value for the width of the Gaussian transverse momentum distribution in pp col- 
lisions. This way our results can be compared to the experimental data on measurable initial 
dijet transverse momentum in a simple way. This analysis has been accomplished [44], using 
recent data on di-hadron correlation [45] . The measured and theoretically extracted momentum 
imbalance value agrees very well with the one applied in our calculations. 
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p + p -> 7T°+X at s 1/2 =200 GeV 

LO PQCD, k=2/3, <k T 2 >=2.5 GeV 2 

NLO PQCD, /t =4/3, <k T 2 >=0.0 GeV 2 

NLO PQCD, k=4/3, <k T 2 >=2.5 GeV 2 
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Fig. 1. Invariant cross section of pion production in p + p — » tv° + X at ^/s = 200 GeV. Upper panel 
displays spectra: LO result at scales Q — (2/3)pt/z c , Qf = (2/3)pr and (fcy) = 2.5 GeV 2 /c 2 (dashed 
line), NLO results at scales Q — Qr — (4/3)pr / z c , Qf = (4/3)pr without intrinsic kr (dash-dotted 
line) and with (fcy) = 2.5 GeV 2 /c 2 (solid line). Lower panel shows the data/pQCD ratios for LO (open 
squares), NLO without (triangles) and with (solid squares) intrinsic kr at the above parameter values. 



4 The parton model for dAu collisions 



Considering the dAu collision, the hard pion production cross section can be written as an 
integral over impact parameter b, where the geometry of the collision is described in the Glauber 
picture: 

E.*£l . /d-M", tAr) M|b - r\) - E, , (4) 

dp J dp 

where the proton-proton cross section on the right hand side represents the cross section from 
eq. (1), but with the broadened widths of the transverse- momentum distributions (3), as a 
consequence of nuclear multiscattering (see eq. (5)). Here £4(6) = f dz pji(b, z) is the nuclear 
thickness function (in terms of the density distribution of the gold nucleus, pau), normalized 
as J d 2 btAu(b) = Aau = 197. For the deuteron, one could use a superposition of a pAu and a 
nAu collision, or a distribution for the nucleons inside the deuteron. For a first orientation we 
follow Ref. [27] in this regard, and apply a hard-sphere approximation for the deuteron with 
A = 2 for estimating the nuclear effects. Also, since ir° production is not sensitive to isospin, 
we continue to use the notation "pA" when talking about the interaction of any nucleon with 
a nucleus. 
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The initial state broadening of the incoming parton distribution function is accounted for 
by an increase in the width of gaussian parton transverse momentum distribution in eq. (3): 

{kj.) pA = (fcr)pp + C ■ h pA (b) . (5) 

Here, (k^) pp is the width of the transverse momentum distribution of partons in pp collisions, 
h p A(b) describes the number of effective nucleon-nucleon (NN) collisions at impact parameter b, 
which impart an average transverse momentum squared C. The effectivity function h p A(b) can 
be written in terms of the number of collisions suffered by the incoming proton in the target 
nucleus, h p A(b) = VA{b) — 1. Here, va(P) = &NNtA{b) , with ctmn being the inelastic nucleon- 
nucleon cross section. Our preliminary results on the determination of the factor of C and v(b) 
in an analysis in NLO confirm the findings of Ref. [11], where the systematic analysis of pA 
reactions was performed in LO and the characteristics of the Cronin effect were determined 
at LO level. Following Ref. [11], we assume that only a limited number of semi- hard collisions 
(with maximum VA{b) ma x = 4) contributes to the broadening, and the factor C = 0.4 GcV 2 /c 2 . 
To give an indication of the dependence of the results on the value of C, we carried out a sample 
calculation with C = 1 GcV 2 /c 2 . 

At RHIC energies, parton broadening can be understood in terms of parton-level collisions. 
Elastic or small-angle-inelastic parton-parton collisions become responsible for the modification 
in the transverse momentum distribution, and this extra broadening will be superimposed to 
the original (k^) pp value [10,46]. The above effectivity function, h p A(b), can be connected to 
the effective partonic collision length. In this framework formation time effects may result in a 
saturation-like behavior. In the following, we apply eq. (5) to describe the transverse momentum 
broadening of partons and use the above parameters. 

It is well-known that the PDFs are modified in the nuclear environment. This is taken 
into account by various shadowing parameterizations [47-50]. In the present work, we display 
results obtained with the EKS parameterization, which has an antishadowing feature [48], and 
with the updated HIJING parameterization [47], which incorporates different quark and gluon 
shadowing, and has an impact-parameter dependent and an impact-parameter independent 
version. The impact-parameter dependence is taken into account by a term oc (1 — b 2 /R\), 
which re- weighs the shadowing effect inside the nucleus. 

The impact parameter dependence of the shadowing function can influence Jjif produc- 
tion in dAu collisions in different rapidity windows [51]. The investigation of the centrality 
dependence of the pion production could contribute to the study of this interesting question. 



5 Results on pion production in dAu collisions 

Including the multiscattering and shadowing effects summarized in the previous section one can 
calculate the invariant cross section for pion production in dAu collision. Moreover, introducing 
the nuclear modification factor RdAu, as 

E^da^/d^p _E 1 ,da d , Au { with nuclear effects)/d 3 p 
dAu ~ Nun ■ E^da^ p /d 3 p ~ E^dui Au { no nuclear eff ects)/d 3 p ' 

nuclear effects can be investigated clearly and efficiently, using a linear scale. The value of Nu n 
can be determined using the Glauber geometrical overlap integral as in eq. (4). Here we apply 
the right-most equation, which docs not require the determination of Nun from the Glauber 
model. 

Figure 2 summarizes our results for the nuclear modification factor RdAu in minimum bias 
dAu collisions at y/s = 200 AGeV. In the top panel the LO case with finite intrinsic fcr is 
displayed. The center panel shows the NLO case without intrinsic kr in the pp collision, corre- 
sponding to a well-focused initial beam of partons, but including the nuclear broadening effect 
connected to multiscattering. The bottom panel shows the NLO results with initial intrinsic fcx 
in the pp collision and nuclear broadening. 
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Fig. 2. The -RdAu nuclear modification factor for minimum bias dAu collision in LO (upper panel), in 
NLO without (middle panel) and with (lower panel) intrinsic for- For details see text and Fig. 1. The 
dotted lines indicate the enhancement in RdAu connected to multiscattering only, which is characterized 
by the average broadening per AW collision, C = 0.4 GeV 2 /c 2 . Thick solid lines show the results after 
additionally including a b- independent shadowing parameterization from the HIJING model [47]. The 
thin line in the bottom panel was obtained with C = 1 GeV 2 /c 2 . The dash-dotted lines display the 
influence of the EKS-shadowing [48], which has an anti-shadowing contribution, increasing RdAu- 



The dashed lines indicate the Cronin enhancement connected to the nuclear multiscattering 
characterized by C = 0.4 GeV 2 /c 2 . The peak structure with an enhancement of 25% is clearly 
seen in the case of large initial intrinsic kx (top and bottom panels) , but in the middle panel 
the Cronin peak is shifted to too small values of pr, out of the range of our pQCD calculations 
(this property was discussed in Ref. [52]). 

The influence of shadowing is indicated by thick solid lines in Fig. 2: a relative decrease of 
RdAu is obtained using a 6-independent parameterization of shadowing applied in the HIJING 
model [47]. Altogether, an upto 10% enhancement can be seen in the 3 GeV/c < pr < 5 
GeV/c region, which disappears at higher px- The height of the peak depends on the value 
of the multiscattering parameter C at a fixed shadowing parameterization. To illustrate this, 
we included a calculation with C = 1 GcV 2 /c 2 in the bottom panel of Fig. 2 (thin line). The 
position of the peak is related to the intrinsic kx value in the pp collision: smaller value leads 
to a shift to smaller px [52]. The EKS parameterization [48], which contains anti-shadowing, 
leads to a small surplus in RdAu, as expected (dash-dotted lines). 

Figure 3 displays the centrality dependence of our NLO results (n = 4/3, (k T ) pp = 2.5 
GeV 2 /c 2 , C — 0.4 GeV 2 /c 2 ) for the nuclear modification factor RdAu with the HIJING shad- 
owing [47]. The left column shows the 6-independent cases. Here, the central bins (upper panel) 
yield essentially overlapping results, because both, multiscattering and shadowing depend on 
the length of the target matter (approximately constant in these cases) . Shadowing finally wins 
moving toward peripheral collisions (lower panel), because the Cronin effect disappears with 
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Fig. 3. The impact parameter dependence of the nuclear modification factor RdAu for b-independent 
shadowing (right column) and 6-dependent shadowing (left column). The lines correspond to different 
centrality bins with an increment of 10%. The 5 "central" curves (0-50%) are shown in the upper 
panels, and the remaining 5 ("peripheral") cases (50-100%) can be seen in the lower panels. 

decreasing u^(b), but a slight shadowing is generated even in the most peripheral cases in the 
recent parameterization of Ref. [47]. However, there is a strong suppression in central collisions 
in the 6-dependent parameterization, which is not balanced by the Cronin effect (see upper right 
panel). Shadowing is rapidly decreasing in this parameterization at larger impact parameter 
(bottom right panel), and the Cronin effect becomes dominant for a certain &-range, increasing 
the nuclear modification factor. Finally, in peripheral collisions the Cronin effect vanishes, and 
in lack of strong shadowing the nuclear modification factor recovers to unity at a reasonable 
transverse momentum value. Precise data on the centrality dependence in the nuclear modi- 
fication factor could yield information about the interplay between multiscattering (including 
saturated vs. non-saturated parameterizations of multiscattering) and shadowing in the dAu 
collision. 

Figure 4 is constructed to emphasize the structure and properties of Fig. 3. The centrality 
dependence of RdAu is shown at fix px = 3.7 GeV/c, in the peak region of RdAu- It is easy to 
appreciate the difference between the 6-independent (solid line) and the fe-dependent shadowing 
(dashed line) from this Figure. 

Figure 5 (top panel) displays our results for the nuclear modification factor in minimum bias 
dAu collision for 6-independent and ^-dependent shadowing parameterizations from HIJ1NG. 
Here we have k = 4/3 for the scales in the NLO calculation, and we use (k^) pp = 2.5 GeV 2 /c 2 
and C = 0.4 GeV 2 /c 2 , as in the bottom panel of Fig. 2. In spite of the very different impact 
parameter dependence of RdAu the minimum bias results are very close to each other. This is 
because many details are averaged out in minimum bias data, and the final result is no longer 
sensitive to the 6-dependence of shadowing. These results can be directly compared to the min- 
imum bias experimental data. In minimum bias data we expect a larger enhancement in 
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Fig. 4. The centrality dependence of the nuclear modification factor RdAu in the peak region, at fix 
pr = 3.7 GeV/c. The solid line corresponds to the case of fe-independent shadowing, the dashed line 
describes 6-dependent shadowing. 



the nuclear modification factor at this energy, due to the presence of protons and antiprotons, 
which are known to yield anomalous p/ir ratios [53]. Uncertainties in the treatment of fragmen- 
tation into protons [54] preclude a detailed prediction of RdAu f° r at the present time. On 
the other hand, parton fragmentation may overlap with parton coalescence in the intermediate 
transverse momentum range, as it was investigated in AuAu collisions [55-57], and this effect 
may have influence in dAu collisions. 

In the bottom panel of Figure 5 we display the double ratio of the nuclear modification 
factors for central to peripheral collision, R* = RdAu(0 — 50%) / RdAu(50 — 90%). The full line 
indicates the ratio R* for the 6-independent shadowing and the dashed line corresponds to the 
6-dependent shadowing. In the window 2 GeV/c < pr < 3 GeV/c the obtained difference is 
~ 30%, which could be seen in case of high precision data are available. 

Figure 6 displays recent high precision experimental data at different centralities on pion 
production in dAu collisions at ^J/vw = 200 GeV [3]. We applied our model and calculated 
pion production, including HIJING shadowing (solid lines) and EKS shadowing (dash-dotted 
lines). Although the size of the error bars does not allow us to draw a firm conclusion about 
centrality dependence, selecting transverse momentum windows at px ~ 3.7 GeV/c, we may 
see a tendency similar to the dashed line in Fig. 4. However, in a recent paper [58], applying 
the HKM shadowing function, we investigated in detail the most central dAu data and have 
found an indication of the presence of jet energy loss in cold matter. This finding means that 
theoretical descriptions have to deal with nuclear shadowing, multiscattering, and jet energy 
loss already in dAu collisions, and may possibly lead to a structure similar to the full line in 
Fig. 4. This interplay can be investigated quantitatively if the precision of the data will become 
even higher by the analysis of a forthcoming dAu run. 
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Fig. 5. Top panel: the nuclear modification factor RdAu for minimum bias dAu collisions obtained from 
our NLO calculation in case of b- independent (solid line) and 6-dependent shadowing (dashed line). 
Experimental data are from the PHENIX Collaboration [24]. Bottom panel: the double ratio of the 
nuclear modification factors for central to peripheral collision, R* = RdAu(0 — 50%) / RdAu(50 — 90%). 
The full line indicates the ratio R* for the ^-independent shadowing, the dashed line corresponds to 
the b-dependent shadowing. 
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6 Discussion 



Summarizing our theoretical results in a NLO pQCD-based parton model calculation for the 
dAu collision, we see a clear enhancement (w 10 — 15%) in the nuclear modification factor 
RdAu after including nuclear multiscattering and HIJING shadowing. Our NLO calculation of 
pion production in the dAu collision confirms earlier LO results, while providing an updated de- 
scription of the invariant cross section and the interplay between multiscattering and shadowing 
effects. We have found that there is a qualitative difference between position-independent and 
position-dependent shadowing; the specific structure of the centrality dependence of the nuclear 
modification factor results in maximum enhancement in semi-central collisions for position- 
dependent shadowing. Thus, the centrality dependence of RdAu may be used to obtain detailed 
information about the effects of multiscattering and shadowing. 

Recent experimental data on tt° production [3] allow us to perform a detailed investigation 
of centrality dependence. Although the minimum bias ir° data display a 10-15% increase in 
the nuclear modification factor supporting the idea of interplay between multiscattering and 
shadowing in the reproduction of the Cronin effect, but detailed conclusion can not be obtained 
at the recent precision of the data. The /i ± data [4,5] appear to have an even larger increase, 
but we do not have a theory comparison for these data at present, because the proton (and 
antiproton) production is not well described by recent pQCD based parton models. 
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Fig. 6. The nuclear modification factor RdAu for pion production in dAu collisions at different cen- 
tralizes measured by PHENIX Collaboration [3] . Results from our NLO calculation are displayed with 
b- independent HIJING shadowing (solid line). Dashed line corresponds to b- independent EKS shadow- 
ing. 
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